Half-ordered state in the anisotropic Haldane-gap antiferromagnet 

Ni(C 5 D 14 N 2 ) 2 N 3 (PF 6 ). 
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Neutron diffraction experiments performed on the Haldane gap materiai Ni(CsDi4N2)2N3(PF6)in 
high magnetic fieids applied at an angle to the principal anisotropy axes reveal two consecutive field- 
induced phase transitions. The low-field phase is the gapped Haldane state, while at high fields the 
system exhibits 3-dimensional long-range Neel order. In a peculiar phase found at intermediate 
fields only half of all the spin chains participate in the long-range ordering, while the other half 
remains disordered and gapped. 
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High-field phase transitions in quantum antiferromag- 
nets (AFs) have recently drawn a great deal of atten- 
tion. Of particular interest is field-induced spin freezing 
exhibited by many quantum spin liquids. The massive 
triplet of low-lying 5 = 1 gap excitations (magnons) 
in such systems is subject to Zeeman splitting by ex- 
ternal magnetic fields. A soft-mode quantum phase 
transition occurs when the gap for one member of the 
triplet is driven to zero. The result is a Bose con- 
densation of magnons. In the presence of magnetic 
anisotropy and weak inter-chain interactions, always 
found in real materials, the magnetized high-field phase 
is a Neel-like state with AF long-range order. The phase 
transition itself and the peculiarities of the high-field 
phase have been studied experimentally in several model 
materials including the Haldane-gap antiferromagnets 
NDMAP (Ni(C5Di4N2)2N 3 (PF 6 )) 1 i 2 i 3 i 4 i 5 i 6 and NDMAZ 
(Ni(C5H 14 N 2 )2N 3 (C104))r^ the bond-alternating 5 = 
1 chain NTENP ( [Ni (iV, JV'-bis ( 3-aminopropyl) propane- 
l,3-diamine(/i-N02)]C104)2i and dimer systems such as 
TlCuClgia and Cs 3 Cr 2 Br 9 ii 

It has been long established that magnetic anisotropy, 
which is particularly important for 5 = 1 materials, has 
a strong impact on the phase transition. The value of 
the critical field H c depends on the relative orientation 
of the applied field and the anisotropy tensor i&liLH A 
very interesting case is that of NDMAP. This compound 
features two equivalent sets of Haldane spin chains with 
non-collinear local anisotropy axes. To date, all exper- 
iments were performed in magnetic fields applied along 
the principal axes of the orthorhombic crystal structure 
that are also the macroscopic magnetic anisotropy axes. 
In these geometries all tilts of the local anisotropy axes 
relative to the field direction are the same for the two 
chain types. In the present study we investigate a less 
symmetric scenario, in which the magnetic field is applied 
in a general direction relative to the crystal axes. We 
find two consecutive field-induced transitions from the 
quantum-disordered spin liquid to ordered Neel phase, 
with a novel "half-ordered" phase in-between. 




FIG. 1: Crystal structure of NDMAP in projection onto the 
(ac) crystallographic plane and the geometry of the present 
experiment. Only the N- (small balls) and Ni-ions (large 
balls) are shown. The local magnetic anisotropy axes x and z 
are tilted in the (ac) plane by a = 16° relative to the c axis. 
The tilt direction is opposite for type-A and type-B chains. 
The field H is applied at an angle ip relative to the c axis 
of the crystal and forms unequal angles <f>\ and 4>b with the 
anisotropy axes z\ and zb, respectively. 



The crystal structure of NDMAP (orthorhombic space 
group Pnmn, a — 18.046 A, b = 8.705 A, and c = 
6.139 A) was described in detail in Ref. 4 . The 5=1 
chains run along the crystallographic c axis. The in- 
chain exchange constant is J = 2.6 meV. Inter-chain 
interactions are much weaker: \J±/J\ < 10~ 3 . Mag- 
netic anisotropy is predominantly of single-ion easy-plane 
type with D/J w 0.25. In addition, there is a weak 
in-plane anisotropy term, and the degeneracy of Hal- 



dane triplet is fully lifted. The gap energies for exci- 
tations polarized along the principal anisotropy axes x, 
y and z are A x = 0.42(3) meV, A y = 0.52(6) meV, and 
A 2 = 1.9(1) meVi^ The anisotropy axes are determined 
by the geometry of the corresponding Ni 2+ coordination 
octahedra and, as mentioned above, do not exactly coin- 
cide with crystallographic directions. Instead, the x and 
z axes of the NiNg octahedra are in the (a, c) crystallo- 
graphic plane, but tilted by a « 16° relative to the a and 
c axes, respectively. There are two types of chains related 
by symmetry, and the corresponding tilt directions are 
opposite. Within each set of chains the Ni 2+ -sites form a 
simple orthorhombic Bravais lattice. On the other hand, 
the two sets of chains are displaced by (0.5, 0.5, 0.5) rela- 
tive to each other. The overall lattice of Ni 2+ ions is thus 
a body-centered one. Due to this geometric frustration, 
the two sets of antiferromagnetic spin chains are magnet- 
ically decoupled at the Mean Field level. The main fea- 
tures of the crystal structure of NDMAP are illustrated 
in Fig.Hl 

Our new neutron diffraction experiments were per- 
formed on a fully deuterated single crystal NDMAP sam- 
ple of approximate mass 100 mg. The data were take on 
the D23 lifting-counter diffractometer at Institut Laue- 
Langevin. Sample environment was a vertical-field cryo- 
magnet with a dilution refrigerator insert. The data were 
taken at T ~ 35 mK in magnetic fields up to 6 T. Un- 
fortunately, there was no possibility to rotate the sample 
in situ during the experiment. Re- mounting the sam- 
ple to explore several orientations was not an option ei- 
ther, as the crystals are known to shatter and deterio- 
rate rapidly during cooling/heating cycles. For this rea- 
son only one experimental geometry was realized, with 
the magnetic field applied at an angle ip — 14.2° to the 
c axis, in the (a, c) crystallographic plane, as shown in 
Fig. ^ F° r one se t °f s pi n chains that we will refer to 
as "type A" , the field was thus almost exactly parallel to 
the main anisotropy axis, the corresponding angle being 
4>s = a — ip = 1.8°. For type-B chains the angle between 
the field direction and the local Ni 2+ anisotropy axes was 
considerably larger, 0b = a + ip = 30.2°. 

Our main experimental result is drawn from the mea- 
sured field dependence of the (0, 0.5, 0.5) and (1, 0.5, 0.5) 
magnetic Bragg intensities at T = 35 mK. These data 
are plotted Fig. [21 In the geometry of the present experi- 
ment two distinct anomalies are detected at = 3.4 T 
and iJ c (2) = 4.1 T. Below hP there is no antiferromag- 
netic Bragg scattering in NDMAP that retains its spin- 
singlet ground state. At antiferromagnetic Bragg 
reflections simultaneously appear at both (0, 0.5, 0.5) and 
(1,0.5,0.5) reciprocal-space positions. These two peaks 
remain of roughly equal intensity in the field range 
H^ 1} < H < Hi 2) . Above the 2nd transition at H ( c 2) the 
intensity of the (0, 0.5, 0.5) peak starts to increase more 
rapidly. In contrast, the (1,0.5,0.5) peak intensity flat- 
tens out and even decreases slightly. Typical scan across 
the two magnetic Bragg reflection collected at H = 6 T 
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FIG. 2: Main panels: field dependence of two magnetic Bragg 
peak intensities measured in NDMAP at T — 35 mK. The 
magnetic field is applied in the (a, c) crystallographic plane 
at an angle 14.2° to the c axis. The solid lines are as described 
in the text. Insets: /i-scans across the two magnetic Bragg 
peaks at H — 6 T. 



are shown in the insets in Fig. [2] 

Though a detailed determination of the magnetic 
structure has not been performed, a survey of several 
magnetic Bragg peaks at H = 6 T revealed a consistent 
intensity pattern. It was found that (ft,+fc+/)-even reflec- 
tions are considerably stronger than the (h + k + l)-odd 
ones. Within each of these reflection classes the inten- 
sity is a smooth function of wave vector transfer, typical 
of the combined effects of the magnetic form factor and 
neutron polarization factors. Such behavior is consistent 
with an almost-collinear antiferromagnetic spin arrange- 
ment on a body-centered lattice, which was previously 
shown to be realized in NDMAP in high magnetic fields 
applied strictly along the high-symmetry crystallographic 
directions 4£ 

The field-dependent behavior observed in the present 
study is in stark contrast that previously seen for H\\c 
(i.e., (f>A = 0b), where only a single transition was de- 
tected. The unusual two-stage transition can be under- 
stood in the framework of a simple model first proposed 
by I. Harada^ The key idea is that long-range ordering 
occurs in each set of spin chains independently 

The critical fields for an individual spin chain in an 
arbitrarily directed magnetic field can be estimated us- 
ing perturbation theory 17,18 Though clearly oversimpli- 
fied, for a magnetic field applied parallel to any of the 
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principal anisotropy axes this method is known to give 
the same values of H c as more sophisticated calculations 
based on the quantum non-linear sigma-modeU^ or map- 
ping to Mayorana fermions, 12 For a field in the (x, z) 
plane applied at an angle <fi to the magnetic easy axis z, 
the perturbative result for H c is£ 

^bHc = > —— 7-5-7 — 5~7- (1) 

y g£A x sin 4> + giA z cos 2 <f> 

In this formula g z and g x are components of the Ni 2+ 
gyromagnetic tensor. Making use of the previously 
measured gyromagnetic ratios^ for the geometry of the 
present experiment one gets H c = 3.8 T and H c — 4.3 T 
for type- A and type-B chains, respectively. While some- 
what larger than the measured values, these two fields 
can be associated with the two observed ordering transi- 
tions at He ancl hP . 

Below He both types of spin chains are in a quantum- 
disordered gapped state. As the external field exceeds 
HP at zero temperature, individual type-A chains ac- 
quire long-range Neel order. Weak interactions be- 
tween type-A chains stabilize this ordered state at non- 
zero temperatures and ensure the existence of true 3- 
dimensional long-range static antiferromagnetic spin cor- 
relations. Nevertheless, considering that inter-chain in- 
teractions are very weak, in the field range HP < 
(2) 

H < He type-B chains must remain in the quantum- 
disordered gapped phase. The corresponding Ni 2+ ions 
carry no static magnetization and do not participate 
in the long-range Neel order. This peculiar phase of 
NDMAP, where half the Haldane spin chains remain 
gapped while the other half participate in long-range an- 
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Fitting this form to the experimental data (solid lines in 
Fig-EJ yields the following parameters: hP = 3.42(1) T, 
HP = 4.07(1) T and — 0.37(2) and mf ] /m^ = 
1.50(4). Overall, the experimentally measured field de- 
pendencies are rather well reproduced by the simple 
model. 

The experimental observations of a two-stage transi- 
tion and half-ordered state in NDMAP highlight several 
important aspects of field-induced spin freezing in highly 
anisotropic gap systems. First, the two-stage transition 
is a dramatic demonstration of the fact that gap energies 



tiferromagnetic order can be described as "half-ordered" 
state. In this regime the only magnetized ions in 
NDMAP are located on type-A chains and form a sim- 
ple Bravais lattice. As a consequence, (0,0.5,0.5) and 
(1, 0.5, 0.5) magnetic Bragg reflections have the same 
structure factor. Their intensities should differ only 
slightly due to slightly different form- and polarization- 
factors. These intensities are proportional to the square 
of the staggered magnetization on the A-sublattice: 

7(0,0.5,0.5) ^(1,0.5,0.5) K I to a| 2 - 

(2) 

The situation changes at He when the gap in type-B 
chains closes as well, and they too acquire static long- 
range AF spin correlations. Now static magnetic mo- 
ments are located on both A- and B-sublattices and form 
a body-centered structure. In spite of geometric frus- 
tration, a definitive relative alignment between spins on 
the two sublattices is established via dipolar interactions 
and/or order- from-disorder fluctuation effects. The two 
magnetic Bragg peaks are then no longer equivalent. As- 
suming a collinear alignment of A- and B-type spins, 
their intensities are given by /(o, 0.5,0. 5) x l m A + ^b| 2 
and /(i,o.5,o.5) l TO A ~ ^b| 2 . As both staggered magne- 
tizations increase with field, /(o, 0.5,0. 5) increases rapidly 
and 7(1,0.5,0.5) levels off. 

To construct a phenomenalogical semi-quantitative de- 
scription we can assume that wia = m^ \H/Hc ~ 
lfe(H-HP) and m B = rog \H/hP - 1\^9(H^hP). 
Assuming that the ordered moments on the two chain 
subsets are collinear (and thus the corresponding po- 
larization factors for neutron scattering intensities are 
equal) , the expression for the measured Bragg intensities 
can be written as: 



(2a) 



(2b) 

I 

and transition fields are influenced by local anisotropy pa- 
rameters, rather than macroscopic magnetic anisotropy 
routinely measured with bulk methods. The second con- 
clusion is that the principle driving force of the phase 
transition is the tendency of individual spin chains to 
form long-range Neel order at T = in sufficiently strong 
fields. Inter-chain interactions are of course needed to 
stabilize this order in 3 dimensions at T > 0, but play 
only a minor role in determining the actual transition 
field at T — > 0. In future experiments it would be very in- 
teresting to investigate the angle-dependence of the tran- 
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sition fields in NDMAP in more detail. 
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